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ABSTRACT 

Massive stars lose mass in the form of stellar winds and outbursts. This material 
accumulates around the star. When the star explodes as a supernova (SN) the resulting 
shock wave expands within this circumstellar medium. The X-ray emission resulting 
from the interaction depends, among other parameters, on the density of this medium, 
and therefore the variation in the X-ray luminosity can be used to study the variation 
in the density structure of the medium. In this paper we explore the X-ray emission 
and lightcurves of all known SNe, in order to study the nature of the medium into 
which they are expanding. In particular we wish to investigate whether young SNe 
are expanding into a steady wind medium, as is most often assumed in the literature. 
We find that in the context of the theoretical arguments that have been generally 
used in the literature, many young SNe, and especially those of Type Iln, which 
are the brightest X-ray luminosity class, do not appear to be expanding into steady 
winds. Some Ilns appear to have very steep X-ray luminosity declines, indicating 
density declines much steeper than r~ 2 . However, other Ilns show a constant or even 
increasing X-ray luminosity over periods of months to years. Many other SNe do not 
appear to have declines consistent with expansion in a steady wind. SNe with lower 
X-ray luminosities appear to be more consistent with steady wind expansion, although 
the numbers are not large enough to make firm statistical comments. The numbers do 
indicate that the expansion and density structure of the circumstellar medium must 
be investigated before assumptions can be made of steady wind expansion. Unless a 
steady wind can be shown, mass-loss rates deduced using this assumption may need 
to be revised. 

Key words: circumstellar matter; stars: massive; stars: mass-loss; supernovae: gen- 
eral; stars: winds, outflows; X-rays: ISM 



1 INTRODUCTION 

In 1972 , in a r eview article on supernova remnants (SNRs), 
IWoltierl (|l972l ) did not devote much attention to the ambi- 
ent medium into which the SN was expanding, suggesting 
that the "importance of the swept-up matter is compara- 
tively minor" in the ejecta-dominated stage, as " everything 
depen ds on the details of the explosive process" . IChevalieri 
1 19771 ) also considered the interaction of SNRs only with the 
interstellar medium. However, by 1982 the idea appears to 
have been established that core-collapse supernovae (SNe) 
expand not into the interstellar medium but the wind-blown 
circumste l lar med ium (CSM) ejected by the progenitor star. 
IChevalieri (|l982al lbh derived self-similar solutions for this in- 
teraction. Since then this concept has been widely tested and 
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accepted, so much so that nowadays it is commonly accepted 
that young core-collapse SNe expand into the wind-blown 
structures created by their progenitor stars, and the result- 
ing X-ray emission is due to the forw ard and reverse shocks 
inter acting with the ambient medium (|Chevalier fc Franssonl 
Il994h . 

The wind from a star is defined basically by two pa- 
rameters, the mass-loss rate and the wind velocity. If these 
parameters are constant it is referred to as a steady wind. 
Another common paradigm that frequently now appears 
in the literature is that young core-collapse SNe are ex- 
panding into steady winds. This has been used to compute 
the mass-loss rates of X-ray supernovae (|lmmler &: Lewinl 
120031: llmmler fc Kuntj 120051). and especiall y optical SNe 



zUUo; llmmler fe Kuntz zUU o]), and es pecially 
jSmith et all 120081 . l2010l : iKiewe et all l20ld ). In very few 
cases has an attempt been made to determine whether the 
wind is in fact steady, and some of the calculations were 



2 Dwarkadas & Gruszko 



shown to be inconsistent i n this respect with the results 
obtained i|Dwarkadasll201ll ). However, a steady wind is a 
reasonable approximation that makes analytical calculations 
possible, whereas the alternative is to do much more compli- 
cated calculations which may not even be possible without 
substantially more data. 

A steady wind with a constant mass-loss rate (M) and 
wind velocity (v w ) results in a density p — M / (Aivr 2 v w ) . 
The density profile thus evolves as r~ 2 . If the wind param- 
eters vary with time, the density profile would deviate from 
an r~ 2 evolution. Thus by determining the nature of the 
density profile into which young core-collapse SNe are ex- 
panding, one could infer the variation in wind parameters 
(if any), and the nature of the progenitor mass- loss. Super- 
nova expansion into the surrounding medium allows a nat- 
ural probe into the density structure of this medium. Since 
the SN expansion velocity can be about 10-1000 times the 
velocity of the surrounding medium, depending on whether 
the progenitor was a Wolf-Rayet (W-R) star or a red super- 
giant (RSG) star, the SN ejecta allows us to probe 10-1000 
years of wind evolution in 1 year of SN evolution. 

The X-ray emission resulting from the circumstellar in- 
teraction is proportional to the density structure (see Sj2}, 
and therefore the evolution of the X-ray light curve can be 
used to probe t he variations in de nsity structure in the am- 
bient medium. IChevalierl (|l982bl) studied the X-ray emis- 
sion from core-collapse SNe, and formulated an expression 
for their X-ray luminosity in terms of the evolutionary pa- 
rameters. This has since been use d to describe the X-ray 
light curves of SNe. ISchlegell (| 19951 ) reviewed X-ray observa- 
tions of all 5 X-ray SNe known till that time. The pace of 
discovery has picked up with the advent of Chan dra, XMM- 
New ton and Swift, with llmmler fc Lewir] (120031) list i ng 15 
SNe,lSchlege]|2006) having 25 on his list, and llmmler] (|2007T ) 
listing 32 SNe.Q 

In this paper we aim to use published X-ray light-curves 
of SNe to explore the nature of the medium into which they 
are evolving. We will explore what can be learned from X- 
ray light curves in 321 an d the limitations of this technique. 
In the next section [3] we show the lightcurves of X-ray SNe 
and fits to the data, and discuss the resulting implications. 
We analyze these results for individual SNe and SNe types 
in 3U Finally, in 39 we summarize our results and discuss 
their consequences. 



2 WHAT CAN WE LEARN FROM X-RAY 
LIGHT-CURVES OF SNE? 

A very thorough discussion of the evolution of the X- 
ray light curv e from young SNe was carried out by 
iFransson et al.1 (|l996l ). in reference to SN 1993J. In this 
discussion, and throughout this paper, we will make use 
of many of the results described therein. For the sake of 
completeness, we present here a simplified analysis with ba- 
sic results that develops and summarizes the most impor- 
tant ideas. This s umm ary follows the description given in 
iDwarkadas et al.l l|2010l ). 

1 An extremely useful compilation of X-ray SNe, maintained 
by Dr. Stefan Immler, is available at the following web site: 
lheawww.gsfc. nasa. gov / users / immler / supernovaejist . html 



The expansion of a SN shock wave into the ambient 
medium gives rise to a forward and reve rse-shocked struc - 
ture, separated by a contact discontinuity (IChevalierlll982al ) . 
The shocks heat up the medium to high temperatures, pro- 
ducing X-ray emission. The X-ray luminosity L x from a 
source depends on the electron density n e , the emitting vol- 
ume V and the cooling function A 

L x ^ n . 2 AV (1) 

The density of the medium into which the SN shock 
wave is expanding is assumed to go as r~ s , with s=2 for a 
steady wind. The emission arises from a thin shell of radius 
Ar in between the forward and reverse shocks, at a mean 
radius r , whose volume V can be expressed as 4-7rr 2 Ar. Al- 
though the wind may not have been steady, the evolution 
will be self-similar as long as both the SN ejecta and the cir- 
cumst ellar medium den sity profiles are described by power- 
laws (jChevalierl Il982ah . In the self-similar case, Ar oc r, 
and therefore V oc r 3 . For a SN in the early stages, the 
post-shock temperature is going to be much larger than 
10 7 K, and the cooling funct ion is assumed to vary as T ' 5 
IjChevalier fc Franssonlll994r ). For a strong shock, T oc v s 2 , 
and therefore the cooling function A oc v a oc r/t in the self- 
similar case. 

Therefore we get 



-2s r 3 
~u r - r 
t 

which gives 



t 



(2) 



(3) 



For s=2 this gives the well-known result that the emis- 
sion decreases inversely with time, L x oc t . This re- 
sult is inherent in the formula for mass-loss derived by 
llmmler fc Lewinl (120031 ). 

In the self-similar case, with p e j oc A t~ 3 v~ n and p cs oc 
r~ 3 , the radius of the SN will evolve as (|Chevalierlll982al ) 

R S Nctt a a = (n - 3)/(n - s) (4) 

Using this in equation [8] gives us that in a general case, 

^ £-(12-7s + 2ns-3n)/(n-s) j-g-j 

This is a general formula that has often been used to 
describe the evolution of the X-ray luminosity. An important 
caveat must be pointed out here - this formula considers the 
total X-ray emission over all temperatures. However X-ray 
satellites measure the emission only in a very narrow band. 
R OSA T covered the region between 0.4 - 2.4 keV. Chandra 
and XMM-Newton cover a slightly larger band, around 0.4- 
10/15 keV, but the effective area is much higher at lower 
energies. Therefore, although equation [5] shows how the X- 
ray luminosity will vary with time over the entire X-ray 
range, it is doubtful that this variation could be observed. 
We observe the variation only for that part of the luminosity 
that falls within a given X-ray mission's observable band. 
This would have the same time dependence if and only if 
the ratio of the luminosity in that discrete band to that of 
the total X-ray region is constant. 

However this is not likely. As the SN expands, we expect 
that the shock will sweep up more material and decelerate. 
The post-shock temperature, proportional to the square of 
the shock velocity, will decrease with time, and the X-ray 
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emission would gradually shift from high to lower temper- 
atures. Since X-ray satellites can detect emission mainly in 
the softer bands, we would expect that the luminosity in the 
lower energy range would increase at the expense of that in 
the higher energy range. Thus, the ratio of X-ray emission 
we observe to the total should increase over time, leading to 
a time dependence flatter than the overall time dependence 
in a steady wind of t~ x . 

In some cases, this can be derived. From above, for a 
non-radiative shock, the shock velocity goes as v 3 h oc t 01-1 , 
and the shock temp erature therefore dec reases as T s h oc 
£ 2 ( a -i) a s shown in lFransson et ail |l996T ). the luminosity 
in a given energy band with E << kT s h will then go as: 



L, — (15 — 6s+ns — 2n) /(n — s) 



(10) 



L X {E) ~ t 



-(6-5s + 2ns-3n)/(n-s) ^ .-0 



(6) 



Compared to equation [5l note that this gives a flat- 
ter time dependence, as expected. For s = 2 we get that 
L x ~ t-(™- 4 )/(™-2). The index /? is now always < 1, im- 
plying that the time dependence we see will not be t _1 but 
will be flatter, and is now a function of the ejecta density 
power-law. For n > 5, as required by the Chevalier solution, 
0.33 < < 1. For commonly used values n=9, (3 = 0.714, 
and for n=ll, /3 = 0.77. Therefore, even for a steady wind, 
we are more likely to see L x ~ £-(0-75±o.05) j t s h omc i De 
noted that for s = 2, we will not see a steeper luminosity 
dependence than t _1 , or flatter than t~ 33 . 

For s / 2 the situation is more complex. For s = 1, 
we get P = — 1, i.e. the flux is linearly increasing with time, 
irrespective of the value of n. Thus a decreasing wind density 
can still give rise to an increasing flux. For s — 1.5 we get 
that /3 = — 1.5/(n — 1.5) is always negative, and therefore 
that the luminosity is increasing with time, whereas the total 
X-ray luminosity from equation[S]is actually decreasing with 
time. Thus this clearly shows the difference between what 
we would observe with X-ray satellites, and what the total 
X-ray luminosity is doing. 

Since we can get both decreasing and increasing lumi- 
nosity dependence with time, it follows that for some spe- 
cific values we could observe a constant luminosity with 
time, i.e. /3 — 0. This depends on the value of the ejecta 
density. For n — 9 we obtain from equation [5] that f3 = 
(13s - 21)/(9 - s). Therefore, for s = 21/13 ~ 1.62 the 
luminosity would be constant with time. 

If the temperature of one of the shocks, especially 
the reverse shock, is lower than about 3 xlO 7 , then 
the cooling function beha vior changes to A oc y -0 - 6 
jChevalier fc Franssonlll994T ). For simplicity and to illustrate 
the difference in evolution, we may assume that the depen- 
dence goes as T -0 ' 5 . Then we get, following the argument 
above, that the A dependence can be inverted, giving 



-2s t 3 
r 

which leads to 



(7) 



L x ~ r*- M t (8) 
and therefore, assuming that r a t a , we get 

L X (E) ~ t -( 2 ^-5 s -3"-+6)/(™-«) ^ t ~P ( 9 ) 

The situation becomes more complicated if the shock 
becomes radiative. This has also been discussed by 
Fransson et al l l| 19961 ). In this case the luminosity goes as 



while, as long as the cool shell is optically thick, the 
luminosity in a band evolves as 



L X (E) ~ y0.16 i (3-2s)(n-3)/(n-s) 



(11) 



where the luminosity is now proportional also to the 
electron temperature T e with time. It is interesting to note 
here that for a density gradient s — 1.5 the flux is con- 
stant irrespective of the value of the ejecta gradient n. For 
s > 1.5 the luminosity is decreasing, while for s < 1.5 it 
is increasing. Note that for a steady wind with s = 2 the 
time dependence /3 = a, i.e. the luminosity evolution has 
the same time dependence as the radius evolution without 
cooling, although with the opposite sign. 

We emphasize that the light curves of young SNe ex- 
panding into a steady wind, observed in a narrow band, 
would not be expected to have the ca nonical t" 1 depen- 
dence that so many have postulate d (Pooley ct al. 20021; 
l lmmler et al.ll2002l; llmmler fc Lewinl^OO^Ilmmler fc Kuntzl 
120051 ; ISoria fc Perna 2008), but would have a natter depen- 



dence that is a function of the ejecta power-law . This has 
been very nicely shown by [Fransson et all i| 19961 ) but is rou- 
tinely ignored. 

Even though this approximate theory works nicely in 
some cases, it is important that we remain cognizant of its 
limitations. This includes the fact that a power-law profile 
is assumed for both the ejecta and the surrounding wind 
medium; that self-similarity is assumed; and that the same 
cooling function is implied for both forward and reverse 
shocks. It is possible that one or both of the ejecta and 
CSM density profiles do not decrease as a power-law; that 
the solution does not depict self-similar behaviour, and that 
different temperatures result in different cooling functions at 
each shock. While in the later case the result for each shock 
could possibly be worked out in a similar manner, the total 
X-ray emission from the SN would depend on details such as 
how much was emanating from each shock in which temper- 
ature range, and would not be easy to compute analytically. 
Finally, such a theory computes the emission due only to cir- 
cumstellar interaction, and does not consider X-ray emission 
due to other mechanisms such as Inverse Compton scatter- 
ing or a central pulsar, or due to other components such as 
radiative shocks in a clumpy medium. Thus, although there 
is a wide range over these ideas are applicable, it is prudent 
to keep the shortcomings in mind. 



3 OBSERVED X-RAY LIGHT CURVES OF 
YOUNG CORE-COLLAPSE SNE 

3.1 X-ray Data: 

In order to study X-ray SNe light curves, we have compiled 
data on SNe that was available in the literature, using pub- 
lished values of luminosity /flux. The time after outburst is 
also taken from the literature, although in cases where dif- 
ferent times of explosion were listed for the same SNe in 
different papers, we have appealed to the optical light curve 
for clarification. Table 1 lists the various papers from which 
the data were compiled. In some cases the data were not 
explicitly listed in a table but had to be read off a fig- 
ure, which introduced further error in the numbers. In a 
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few cases (SN 2005kd, SN 2006jd), we have downloaded the 
publicly archived data-set and computed the flux, in order 
to procure at least one more data-point so that we may 
construct a light curve. In one case (SN 1993J) the most 
current data point is from a Chandra observation on which 
the first author is PI. We have chosen not to include GRB- 
related SNe in this compilation, as their X-ray emission is 
not expected to be due purely to circumstellar interaction. 
We have also not included SN 2008D, a very well observed 
SN whose early lightcurve at le ast was attributed to shock 
breakout (ISoderberg et al]|2008l ) and thus did not fit in with 
the circumstellar model discussed herein. 

Figure[T]plots the datapoints for all X-ray observed SNe 
that we found in the literature with at least 1 published data 
point. The number is only 42, although it is rapidly increas- 
ing, thanks mainly to Swift. We plot light curves for almost 
all published SNe with multiple data points, enabling an 
X-ray light curve to be constructed. As far as possible we 
have constructed the light curve in the 0.4-2.4 keV or re- 
lated energy range, which usually has the most data for SNe 
that exploded during or before the ROSAT era, and thus 
ensures the largest duration light curves. In some cases, as 
mentioned in Table 1, a larger energy band is used. Most of 
these are SNe observed with Swift. All points on a lightcurve 
lie in the same energy range, which is given in Table 1. In 
order to accomplish this using only published data, we had 
to select the energy range in which the maximum number 
of SN fluxes were quoted, which sometimes meant not using 
available flux results lying outside the energy range. Given 
our intent of assessing what the published light curves are 
telling us, no attempt was made to interpolate or extrapo- 
late published data from a given energy range to another, 
which would have required assuming a spectral model that 
was generally not available. 

While the number of SNe is gradually increasing, it is 
still not very large, with only 20 SNe found in the literature 
which have flux values at more than one epoch (although 
large amounts of unpublished data exists in the archives). 
Given the current availability of Chandra, XMM-Newton 
and especially Swift, we hope that this situation can be rec- 
tified in future. Only one SN (1993J) has been observed 
regularly since its birth, and has the most complete light 
curve, although gaps still exist. SN 1987A is probably the 
most observed X-ray SN, due to its closeness and unique 
nature. Hard X-ray data exist in the first 1000 days from 
Ginga (|lnoue et al.l fl991). with progressively better resolu- 
tion data thereafter from all the major X-ray observatories, 
and almost every X-ray telescope, including high-resolution 
gratings. 

The plot shows that very few SNe possess a monotonic 
light curve in the given X-ray band. This may be partially 
due to the published error bars being too optimistic, but in 
some cases it is evident that the light curve tends to jump 
around. 



3.2 Fits to the Light Curves: 

Following the discussion in j|2] we fitted the light curves 
of all SNe with multiple data points (Figure [TJ using an 
expression of the form L x = -P(O ) x t p( - x \ The fitt ing was 
done using the MPFIT routine <|Markwardtl 120091 1 in the 
IDL programming language. In Figure [2] we show the fits for 



various SNe. The legend gives the value of the parameter 
P(l) in the expression above, which corresponds to — /3 from 

m 

In many cases the use of a simple power-law fitting func- 
tion does not seem appropriate, as is evident not only from 
the large values of x-squared that were found, but merely by 
comparing the fit to the data points. It appears that a more 
complicated function is needed, or perhaps a break is needed 
in the power-law. In one case, SN 1978K, the fit seems en- 
tirely off, and it is clear that something more sophisticated 
is needed. 



3.3 X-ray lightcurves by SN Type: 

In order to study the variation of the X-ray emission with 
type of SNe, we show in figure [3] the light curves of all SNe 
plotted as a function of the type of SN. It is clear from 
this plot that Type Iln SNe show the highest luminosities, 
and Type IIP SNe the lowest. These bracket the Type IIL's, 
Type lib's and Type Ib/c. We caution that we have not 
included GRB-SNe here, and while the sample is probably 
enough to make general inferences, it is still not large enough 
to draw detailed statistical conclusions. 



4 ANALYSIS OF SN X-RAY LIGHTCURVES 

A survey of the fits to the data in figure [2] reveals the di- 
versity in the X-ray expansion of SNe. About half the SNe 
decay with f3 < 1, a necessary (but certainly not suffi- 
cient) condition if they are expanding in a steady wind. 
Two of them (SN 1987A and SN 1996cr) show f3 < 0, 
displaying X-ray emission that appears to increase over 
long time-scales, as has been well documented in the past 
Chevalier fc Dwarkadasl Il995l; iMcCrav) Il993t iPark et al ' 



2006; IMcCrav! 120071 ; iBauer et all 120081 ; iDwarkadas et al 



20101 ). Another two (SN 1979c and SN 2006jd) appear to be 



relatively constant, although we caution that for SN 2006jd 
we have only two data points, and are waiting for more data 
to be published. A few SNe have /3 > 1. In the framework 
of the theory outlined above, these could not arise from SN 
expansion into a steady wind. This group consists of 4 Type 
Iln SNe, which tend to show some of the steepest decreases 
in X-ray emission. Others in this category include SN 2003bg 
and 2008ax, which do not have enough time-sampling over 
even a decade to effective judge their luminosity decrease. 
The steepest light curve belongs to SN 1970G, which drops 
as t~ 3,08 . Interestingly, many of these SNe which have steep 
drop-off in light curves appear to be those from whom X-ray 
emission was found only late in their lifetimes, after about 
6 years. In this context, we point out that if the X-ray light 
curve of SN 1993J after 6 years is considered, the fall-off 
would be much steeper than the gent ler overall slope, as has 
been shown in lChandra et all |2009t) . 

We discuss some of our results in more details, in the 
context of either individual SNe or groups of SNe. 



4.1 Type Iln Supernovae: 

In 19 90, a new class of SNe was introduced by ISchlegell 
(1990). The 'n' denotes narrow lines, which distinguish Type 
Iln from other SNe. Many of these SNe (but not all) have a 
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Table 1. The list of core-collapse supernovae for which X-ray data were available and used in 
this paper. The source of the data is listed, along with the energy range for the luminosity. DG11 
refers to this paper. 



SN Name 


Energy Range (kev) 


1986J 


.5-2.5 


1988Z 


.2-2.0 


1995N 


.1-2.4 


1998S 


.4-2.0 


1978K 


.5-2.0 


1993J 


.5-2.4 


1979C 


.3-2.0 


2002ap 


.3-2.0 


1994W 


.1-2.4 


1996cr 


.5-2.0 


1981K 


.5-2.0 


1987A 


.5-2.0 


2010F 


.2-10 


2009mk 


.2-10 


2009gj 


.2-10 


2009dd 


.2-10 


2008ij 


.2-10 


2008ax 


.2-10 


2008M 


.2-10 


2006jd 


.2-10 


2006jc 


.2-10 


2005kd 


.2-10 


2005ip 


.2-10 


2004et 


.4-8 


2004dj 


.5-8 


2002hi 


.5-10 


2002hh 


.4-8 


2001ig 


.2-10 


2001gd 


.3-5 


1999em 


.5-8 


2003L 


.5-5 


2003bg 


.3-10 


1999gi 


.3-10 


1970G 


.3-2 


2006bp 


.2-10 


1941C 


.3-8 


1959D 


.3-8 


1968D 


.3-8 


1980K 


.5-4 


2010jl 


0.2-10 


2010jr 


0.2-10 



Publication 



Houck 
Schlcgcl & Pctrc 
Chandra et al. 
Pooley_e^aL 

Schlegel et al. (2004) ; Lenz fc Schlegel 
Chandra et al. (2009), 
Patnaudc ct al. 
Sutariaetal. 
^chlcgcl 
Bauer ct al. 
Immlcr ct al. 
Hasingcr ct al. (1996); Park et al. 

Russelletal. 
Russell & Immlcr 
Immlcr & Russell 
hnmteretaL 
Immlcr ct al. 
Immlcr 
Immlcr 
Immler et al. (2007) . 
hrrrnleretal. 

Immler et al. (2007) ; Poolev et al. (2007) . 

J^rrmler_fc_Poolcy 

Rho ct al. 
Poole v fe Lewin 



Immler fc Brown 



Poolev & Lcwin 
Poolcy & Lcwin 
Schlcgcl & Ryder 
Perez- Torres ct al. 
Poolev ct al. 
Sodcrbcrg ct al. 
Sodcrbcrg ct al. 

Schlcgcl 
Immlcr fc Kuntz 
(2006) ; Immler et al. 
Scria&Pema 
Soria fc Pcrna 



Soria fc Pcrna 



Schlegel (1995 



Soria fc Pcrna 
Immlcr ct al. 
Immlcr ct al. 



(2005) 
(2006) 
(2005) 
(2002) 
(2007) 
DG11 
(2011) 
(2003) 
(1999) 
(2008) 
(2007) 
(2007) 
(2010) 
(2010) 
(2009) 
(2009) 
(2009) 
(2008) 
(2010) 
DG11 
(2008) 
DG11 
(2007) 
(2007) 
(2004) 
(2003) 
(2002) 
(2002) 
(2005) 
(2002) 
(2005) 
(2006) 
(2001) 
(2005) 
(2007) 
(2008) 
(2008) 
(2008) 
(2008) 
(2010) 
(2010) 



high radio and X-ray luminosity. In fact, as is evident from 
Figure [3] as a class they have the highest X-ray luminosity 
of all SN types. 

As shown in Figured 4 of the Type Ilns (SN 1986 J, SN 
1988Z, SN 1998S and SN 1995N) have a very steep luminos- 
ity decline. SN 1986J X-ray flux decreases as i -1 ' 7 . We note 
here that for 1986J, w e have adopted the flu x values com- 
puted by iHouckl (|2005l ) rather than those by iTemple et al.l 
(|2005l ). We recomputed the values for the Chandra observa- 
tions and find much closer agreement with the former than 
the latter. 

The steep value of the X-ray decline would suggest that 
the SN is expanding in a medium with a density that de- 
clines faster than r~ 2 . SN 1986 J is also well studied i n the 
radio range. VLBI observations l|Bietenholz et al1l2010l ) sug- 
gest that the shock radius is expanding as R sh oc t ' 69 . The 
luminosity decline and radial expansion then provide two 



equations which can be used to determine the slope of the 
ejecta n and the slope of the CSM density profile s. It is 
easy to show however, that there is no finite solution for n 
and s using equations 0] and [5] In other words, there is no 
compatible solution for the case of an adiabatic shock in SN 
1986J. 

If we consider a radiative reverse shock, then it is pos- 
sible to get a finite solution which gives s = 2.72. This ap- 
pears to be consistent with the luminosity decline, which 
indicates a very rapid slope. However, it gives a value of 
n — 3.61, which is no t allowed within the Chevalier solu- 
tion (|Chevalierl [l982al ). Such a low value of n is inconsis- 
tent wit h most computed models of SN outer ejecta densit y 
profiles (|Chevalier fc Sokerlll989l ; iMatzner fc McKeelll999l ). 
That does not completely exclude this solution. The emis- 
sion could arise from the reverse shocked material colliding 
with the flatter part of the ejecta density profile (i.e. the 
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Figure 1. The X-ray lightcurves of all observed X-ray SNe found in the literature, with one or more data points. The observed range 
of luminosities spans about 5 orders of magnitude, from 10 37 to 10 42 ergs s _1 , excluding SN 1987A, which is a very low luminosity SN 
that was observed only because of its proximity. 



self-similar solution is no longer valid). It is also possible 
to envision scenarios such as the shock becoming suddenly 
radiative, implying a transition stage, or expanding into a 
medium that drops steeply in density. 

These results seem to point to a complicated e mission 
struct ure. An alternative solution, as suggested by IChugall 
1 19931 ). is that there are two components to the emission. 
The basic idea is that of a SN expanding in a clumpy wind, 
with the clumps having a low volume filling factor. The ex- 
pansion of the SN in the interclump wind medium is mainly 
responsible for the radio emission and the radial expan- 
sion measured with VLBI, whereas the X-ray emission arises 
from radiative shocks within the clumps. However, this still 
implies that the density of the clumps, or the number of 
clumps, must decrease faster than r~ 2 in order to get the 
appropriate decline. 

Even if the SN is expanding into a clumpy wind, it is 
unlikely that it would have been doing so since it exploded. If 
we assume the current luminosity decline started right from 
the SN explosion, and extrapolate the soft X-ray luminosity 
to 10 days after the explosion, we get that it would be about 
1.5 xlO 44 ergs s _1 . This is an extremely large luminosity, 
especially just in the soft X-ray band, given that we would 
expect the hard X-ray luminosity to be much larger at this 
stage (indeed, at the last Chandra observation, the soft and 
hard luminosities were comparable). Note that if the SN 
really was emitting in X-rays with this luminosity, it would 



have lost all its kinetic energy in a year at that rate. Even 
given the rate of luminosity decline it would have lost all its 
energy in a few years. It is clear that such a high luminosity 
would be unsustainable. 

The first X-ray detection of SN 1986J was made 8.43 
years after explosion, when the luminosity was almost 10 40 
ergs s _1 . If we assume that the SN maintained this con- 
stant luminosity since explosion, it would mean that in the 
first 8.43 years it would have lost a total of 2.5 xlO 48 ergs 
just in the soft X-ray band. The amount lost in hard X-rays 
would most likely be greater since the temperature would 
be higher. The flux at infra-red wavelengths can be up to a 
factor of 100 l arger than that in X-rays, especially at higher 
temperatures l|Dwek et alJl200Sft . thus leading to an improb- 
ably high broadband flux. However it is possible that for SNe 
expanding in a very dense medium, some of the dust may be 
destroyed, and the IR flux may be only a few times larger 
than X-rays. The broad-band flux would still be high, but 
not impossibly so. 

A plausible inference is that the average luminosity of 
SN 1986J over the first 8.43 years was lower than at 8.43 
years, and perhaps orders of magnitude lower. At some point 
before 8 years the X-ray luminosity increased. In terms of the 
clumpy wind model, this could mean that there were fewer 
to no clumps in the first few years, and that the number of 
clumps increased after a few years. An alternative is that 
the SN shock was evolving in a much lower density region 
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Figure 2. Fits to the lightcurves of all X-ray SNe with multiple data points. The slope of the Luminosity v/s time curves, plotted on a 
log-log scale, is given next to each SN. The fits were calculated using the MPFIT routine in IDL (see text). 



earlier, and started interacting with a higher density region 
at a later time. This would imply a sudden deceleration in 
the shock velocity when it began expanding into the high 
density region, which is not seen in the radio VLBI data. 
We will explore the X-ray emission from SN 1986J in detail 
in a future paper. 

SN 1988Z has an X-ray luminosity that decreases even 
faster than that of SN 1986J. Following the arguments 
above, a plausible model appears to be that of a SN shock 
wave interacting w i th a clumpy wind, as suggested by 
IChugai fc Danzigerl (|l994h . Not surprisingly, extrapolating 
back to 10 days, we get an X-ray luminosity of 6.5 x 10 46 
ergs s _1 . At this luminosity, the SN would lose all its energy 
in one day, so its clear that the soft X-ray lightcurve could 
not have had the recorded decline right from the time of ex- 
plosion. Following the line of reasoning above for SN 1986J, 
it appears most likely that SN 1988Z also had a lower X- 
ray luminosity in the early years, which increased at a later 
time. This then further suggests a lower density medium, or 
lack of clumps in the ambient medium, during the first few 
years. 

SN 1995N and SN 1998S also have luminosity declin- 
ing faster than t _1 , although not as fast as 1986J or 1988Z. 
Nonetheless, it does suggest a medium that decreases faster 



than 



Extrapolating SN 1995N luminosity back to 10 



days gives a luminosity of ~ 2.5 x 10 ergs s~ , still larger 
than the observed X-ray luminosity of any known X-ray SN, 
and unlikely (although not impossible), following the argu- 



ments given above for SN 1986J. SN 1998S extrapolated 
back gives a luminosity of ~ 5.8 x 10 41 ergss _1 , approaching 
the upper limit of known X-ray luminosities. In the case of 
1998S the decline is not so steep, and it is possible that the 
luminosity could have declined with this slope since the time 
of explosion. 

The above slopes suggest that the density into which 
the SN is expanding, whether it be due to clumps or a wind, 
falls faster than r~ 2 . While most Type Ilns show a very high 
luminosity, not all of them decline so rapidly. SN 2005kd 
appears to rise and fall, over a time-period of only a few 
100 days. SN 2006jd appears to be more or less constant 
over 2 data points separated by a few 100 days. Neither of 
these would be consistent with an r -2 density in the am- 
bien t medium. SN 1996cr, which has been studied in de- 
tail (IBauer et al. I l2008l ; iDwarkadas et al.ll20ld : iDewev et al.l 



1 |B 

l201ll ) . was increasing in X-ray luminosity for several years, 
before turning over, and is thought to have expanded in a 
low density region surrounded by a dense shell. 

It is clear that although there is great diversity, none of 
the Type Iln appear to be evolving into a steady wind with 
constant parameters. 



4.2 SN 1993J: 

This SN has more X-ray data than any other except SN 
1987A. The X-ray and radio emission has been studied 
in detail by several authors, but no consensus reached 
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Figure 3. The light curves of X-ray SNe plotted as a function of the type of SN. Type Iln SNe are the brightest as a group, therefore 
suggesting that they interact with regions of higher density. Others show much more diversity, with HP's generally having the lowest 
luminosity as a class, although exceptions exist in all classes. 



on what kind of medium the SN is expanding in. One 
of the earliest papers wh ich studied the X-ray light-curve 
|Suzuki fc Nomotol |l995) suggested that the density im- 
mediately around the SN was decreasing as r -1,7 . A sim- 
ilar n3sjdt^_of_j5^ensity slope flatter than r -2 , was found 
by iFransson et~aH (| 19961 ). Following up on a seminal pa- 
per that interpreted the ra dio emission of 1993J in th e 
context of a steady wind (|Fransson fc Biornssonl Il99ct ). 
IFransson fe Biornsso n (2005) asserted that both the X-ray 
and radio lightcurves could be fit by a self-similar expansion 
for SN 1993J, with a power-law ejecta expanding into an 
r~ 2 medium. They critic ized the hydrodynamic modelling of 
ISuzuki fc Nomotol (|l995h for using an ejecta profile that wa s 
hydrodynamically unstable. However. I Nvmark et alj (2009), 
citing the same hydrodynamic models, suggested that the 
X-ray emission was coming partly from an adiabatic shock 
expanding into CNO-enriched ejecta (which was not de- 
scribed simply by a power-law), and partly from a radiative 
shock, with no cont ribution from the circ umstellar shock. 
The calculations of IChandra et al.l (j2009h were consistent 
with the fact that the emission was dominated by a reverse 
shock component, but did not need two emitting compo- 
nents, although they did not study the spectra. The radio 
and VLBI data has always been c ontentious, and th e two 
different views are summarized in iBartel et al] l|2007l ) and 
iMartf-Vidal et al l (|201ll ). 

The soft X-ray light-curve for the first 700 days is well 



fit by a slope @ = 0.34. This could fit an s = 2 solution if 
n = 5, but s uch an ejecta profile is in consistent with explo- 
sion models l|Suzuki fc Nomotol Il995l ). Otherwise it cannot 
be fitt ed with a steady wind expansion. As IFransson et al.l 
(|l996h have shown, for s ~ 1.7 the fit is good with a radia- 
tive shock solution. However, the behaviour after 700 days 
is much more complex and does not appear to fit a constant 
value of ft, thus suggesting that the slope is c hanging at this 
time. A changing density slope was used in IChandra et ahl 
(2009) to fit the hard X-ray light curve. Overall, SN 1993J 
does not appear to fit the expansion pattern of a shock ex- 
panding in a steady wind, especially since the slope of the 
X-ray decay steepens at later times. 

4.3 SN 1979C: 

This SN appears to be expandin g with a more or less c on- 
stant luminosity for several vears. lPatnaude et al.l (|201ll ) in- 
terpret this steady luminosity as evidence for emission from 
a central black-hole component. Although without a detailed 
inspection of the spectra, we cannot validate or dispute their 
assertion, we do wish to mention here that having a more 
or less constant lightcurve is not inconsistent with the cir- 
cumstellar interaction model. In section [2] we have shown 
that for adiabatic shocks for example, n=9, s=1.62 gives a 
constant luminosity in the soft X-ray band, while for radia- 
tive shocks, s — 1.5 gives a constant flux. Therefore having 
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a constant flux over an extended period in SN 1979C (and 
perhaps SN 2006jd) does not contradict the circumstellar in- 
teraction model, although it does indicate that the evolution 
is not into a steady wind. 



4.4 Other SNe: 

SN 2006jc is an unusual SN that experien ced a mass ejec- 
tion about two years be fore becoming a SN (iPastorello et al] 
20071 ; iFolev et alj|2007h . The X-ray emission jlmmler et ail 
20081 ) is consistent with the SN ejecta rising as it impacts 



the dense shell resulting from this mass-loss, followed by a 
decline as the sho ck crosses the dense region (see for instance 
|Pwarkadasll2005h . If it was interacting with an r 2 wind, it 
was only for an extremely short period, but for the most 
part the profile deviates significantly from a steady wind. 

SN 1987A has been very well studied in the litera- 
ture. Its proximity, and the availability of the Hubble Space 
Telescope, has allowed us for once to optically observe 
the interaction. Analysis of the optical, x-ray and radio 
data have shown that the SN evolved into a low density 
medi um before impacting a higher density ionized HII re- 
gion (IChevalier fc Dwarkadas!ll995l ; iPark et al.ll2005l . l200rj . 
l2007l ; iDwarkadasI l2007al lbt) followed by a dense shell, part 



of a bipolar circumstellar nebula. The evolution, and X-ray 
emiss ion, is summarized in reviews bv lMcCravl l|l993l . 12003, 
12007ft . 

SN 1978K also seems to be continually bright in X-rays. 
Its luminosity could be considered constant on average, but 
appears to fluctuate b y factors of a fe w, making it difficult to 
get a good linear fit. Schlegell (|2006h asserts that although 
there may not be too much luminosity evolution, spectral 
evolution does exist. The current X-ray light curve though 
appears inconsistent with a steady wind. 

SN 1970G was discovered only several thousand days 
after explosion. Although it currently has the steepest lu- 
minosity slope /3 = 3.08, it is not clear what the slope was 
over the first 20+ years. Currently the slope does not appear 
con sistent with inte raction with an r -2 medium. 

ISchlegell lj200lT l analysed the X -ray emission from Type 
IIP SNe, and suggested that they would not be strong X-ray 
sources. This is consistent with the lightcurves in figure O 
which shows that Type IIPs have a low X-ray luminosity. 
The fits to the slopes of many of these, as well as other SNe 
not explicitly mentioned above, are theoretically consistent 
with expansion within a steady wind. The closest to that of 
a steady wind is that of SN 19941. In most of these cases, 
there are either too few datapoints, or the fits are not sta- 
tistically good, so it is hard to judge if they are evolving 
in a steady wind ove r an extended time period. However, 
IChevalier et al.l (|2006h have studied the X-ray emission from 
Type HP's assuming they evolve in a steady wind, and their 
calculations tend to match the observations. It must be kept 
in mind that HP's evolve from RSGs, which do not have a 
very high wind velocity, and thus the wind will not extend 
very far out from the star, perhaps a few parsecs at most. 
Thus we would expect that even if the SN shocks were in- 
teracting with winds, at some point early in their evolution 
(few hundred years or less) they would run out of wind ma- 
terial to expand in. 



5 DISCUSSION 

In this paper we have plotted the X-ray lightcurves of al- 
most all published X-ray SNe, in the narrow energy ranges 
in which they were observed. Furthermore, we have fitted 
these lightcurves to a function that goes as t p ^ , and derived 
the value of -P(l), which we have compared to theoretical ex- 
pectations. In many of the cases, and especially for Type Iln 
SNe, we find that the light-curve data in this simple theory 
are not consistent with evolution into a steady wind whose 
density decreases as r ~ 2 . There are of cour s e seve ral limi- 
tations to this theory. IChevalier fc Franssonl l|2006h suggest 
that the x-ray luminosities of type Ib/c SNe can only be 
explained by a non-thermal mechanism, either X-ray syn- 
chrotron or inverse Compton. Our sample does not include 
any GRB related lb/c, and only a few non-GRB lb/c SNe 
that have published data. It is thus possible for scenarios to 
exist where the light curve slopes, although not within the 
range outlined in Sj2l may be consistent with a steady wind. 
Although such exceptions are possible, it would be unusual 
to assume that they would apply to all, or even most, SNe. 
Furthermore, if the latter is indeed true, then it suggests 
that the many tens of papers in the literature that have 
used this theory need to be re-evaluated. 

Our plott ed lightcurves are consistent with those of 
ISchlegell (2006) in those cases w here ove r lap is possible. They 
are inconsistent with those of llmmlerl l|2007l ). who plotted 
light curves of all SNe known till then with just the first and 
last data point, while asserting that "all other X-ray data 
points are along the extrapolated lines". As is evident from 
figure [TJ this does n ot seem to represent th e available data. 
Using that assertion, iKoss fc Immlerl (|2007T ) went a step fur- 
ther, calculating the density profile of the medium into which 
the SN is evolving, and showing that the results are consis- 
tent with a density decline corresponding to a steady wind. 
Inherent in their calculation is that the SN luminosity is de- 
creasing as t , which they take to mean that the ambient 
density is going as r -2 , although as emphasized in Sethis 
is not necessarily true. Therefore, it is not surprising that 
their final result matches that initial assumption. Given the 
arguments above, and the light curves plotted in Figure Q] 
and fitted in Figure [2] we do not agree with their results. 

Often in the literature, the X-ray lightcurve, or even 
a single X-ray datapoint, has been measured, and used to 
compute a mass-l oss rate, assuming t hat the light curve 
decreases as t~ x dPqolev et al.l 12002]; llmmler et al.l 2002 ; 
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Immler fc Lewinlliobll : llmmler fc Kunt 
IKoss fc Immlerl 120071 : ISoria fc Pernal 
2010). As shown above, and most clearly in lFransson et al 
(1996), the f _1 dependence for a steady wind is only valid 
over the total X-ray range, not the narrow bands in which 
SNe are generally observed. Even if it were valid in the nar- 
row range, none of the published lightcurves are actually 
decreasing in luminosity as t , and very few are even close. 
We caution therefore that the resulting mass-loss rates ob- 
tained in this manner from the X-ray lightcurves may be 
incorrect, and not reflective of the true nature of the CSM. 
If, as the lightcurves appear to indicate, the SNe are not 
evolving in an r -2 medium, then the calculated mass-loss 
rates would need to be re-calibrated. Unfortunately, if the 
wind is not steady, then the generally used expressions to 
measure the mass-loss rates are no longer valid, and there 



10 Dwarkadas & Gruszko 



is no simple way to compute the mass- loss rate. Most tech- 
niques usually yield the combination M/v w - disentangling 
the magnitude, and time evolution, of either or both param- 
eters is a difficult task. If the deviation from a steady wind 
is not large, a mass-loss rate obtained by assuming a steady 
wind may be a good approximation early on, but the dif- 
ference increases gradually with increasing radius from the 
center of the explosion. 

An r~ 2 density decline is often assumed in calculations 
of mass-loss rates from opt ic al observations ifemith et al.l 
l200Sl . l201Cl : rKiewe et alJiblch . lKiewe et all (|2010l ) have used 
this in fact to compute the mass-loss rates for 15 type Iln 
SNe, and suggest that they are consistent with LBV pro- 
genitors. It h as not been show n that such a density decline 
is warranted. ISalama nca ( 2003) pointed out many years ago 
that the medium around Type Iln SNe does not necessar- 
ily decrease as r~ 2 ; this asserti o n wa s repeated for specific 
individual cases by [Dwarkadas (201l|). If the nature of the 
density distribution plays such an important part in deter- 
mining the mass-loss rate and thereby categorizing the SN 
progenitor, it is important that the assumption of a steady 
wind density decreasing as r -2 not be made without accom- 
panying proof. 

There exist suggestions from other sources that the cir- 
cumstellar medium in the immediate vicinity of massive 
stars does not have an r -2 density profile, or that such a 
profile doe s not extend f ar out . In their analysis of GRB 
afterglows, ISchulze et al.l (j201lT > find that only a quarter of 
them appear to be evolving in a freely expanding wind, with 
the rest appearing to expand in a constant density medium. 
Furthermore, in all cases except two, they place limits on 
the freely expanding wind region of less than 1 pc, and 
< O.lpc for about half the observed GRBs. Radio obser- 
vations of SNe appear to imply density modulations in the 
circumstellar medium ([Ryder et all |2004| ; ISoderberg et al.l 
l2006l ; IWellons fc Soderberdl201ll i7although its not clear how 
large a deviation from a standard r~ 2 profile is indicated. 

X-ray emission from young SNe, which is mainly ther- 
mal, has the simplest and most direct dependence on the 
density structure of the surrounding medium, and can be 
used therefore to most easily infer this density structure. In 
order to improve this analysis and make more detailed cal- 
culations from SN X-ray lightcurves, it is imperative that we 
have detailed lightcurves, with much better time sampling, 
and over as wide an energy range as possible. Given the 
availability of Chandra, XMM-Newton , Suzaku, and Swift, 
and that they should continue operating for many more 
years, it is urgent that we observe as many SNe as possi- 
ble, and for as long as possible, in X-rays. 
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